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ABSTRACT 

Torsional Alfven waves are likely excited with bursts in rotating magnetars. These 
waves are probably propagated through corotating atmospheres toward a vacuum 
exterior. We have studied the physical effects of the azimuthal wave number and 
the characteristic height of the plasma medium on wave transmission. In this work, 
explicit calculations were carried out based on the three-layered cylindrical model. 
We found that the coupling strength between the internal shear and the external 
Alfven modes is drastically enhanced, when resonance occurs in the corotating plasma 
cavity. The spatial structure of the electromagnetic fields in the resonance cavity is 
also investigated when Alfven waves exhibit resonance. 
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1 INTRODUCTION 

Soft gamma-ray repeaters (SGRs) and anomalous X-ray pul- 
sars (AXPs) have been known as strongly magnetized neu- 
tron stars, magnetars. The hallmark of these objects is to 
repeat X-ray or gamma-ray emissions irregularly. So far, 
four or five known active objects which show frequent X- 
ray emission of tremendously energetic and shorter initial 
bursts (typically E < 10 41 ergs and At ~ 0.5 s) have been 
identified with SGRs (Thompson et al. 2001). Some of them 
occasionally exhibit more energetic events. The giant flare 
1900+14, now associated with SGR 0525-66, was observed in 
1979 (Mazets, Golenetskii, and Gur'yan 1979) for the first 
time and became active again in 1992 (Kouveliotou et al. 
1993) and also in 1998 (Kouveliotou et al. 1998 ; Hurley et 
al. 1999). Recently, the intense 7-ray flare from SGR 1806- 
20 occured on December 27, 2004 was reported (Palmer et 
al. 2005; Hurley et al. 2005; Mereghetti et al. 2005). Rea et 
al. (2005) have investigated the pulse profile and flare spec- 
trum of SGR 1806-20. Their study may potentially give the 
information about the global field structure in the magneto- 
sphere and may further promote theoretical magnetar mod- 
els. Their extream peak luminosity extends up to 10 6 I/Edd 
estimated at a distance of 10 kpc (Mazets et al. 1999; Fe- 
roci et al. 2001). Typical X-ray luminosities of SGRs have 
been measured to be L x = 10 34 -10 36 erg s _1 (Hurley et al. 
2000; Thompson et al. 2000), except with giant bursts SGRs 
1900+14 and 1806-20. Shorter durations of initial flares of 
SGRs are comparable to the Alfven crossing time of the 
core. The energy distribution shows good agreement with 
Gutenberg-Lichter law, which may indicate statistical simi- 
larity to earthquakes or solar flares (Cheng et al. 1996; Gogus 



et al. 1999, 2000). The SGRs have spin periods in a small 
range of P = 5-8 s with rapid spin down rate P ~ 10 -10 
ss _1 , which therefore give characteristic ages P/P ~ 10 3 yr 
(Mazets et al. 1979; Kouveliotou et al. 1998; Hurley et al. 
1999a). 

On the other hand, the physical nature of AXPs is still 
uncertain due to poor observations, but there seem to be 
some similarities and differences between SGRs and AXPs. 
AXPs are energetic sources of pulsed X-ray emission, whose 
periods lie in a narrow range P = 6-12 s, characteristic 
ages P/P = 3 x 10 3 ~4x 10 5 yr and X-ray luminosities 
L x = 5 x 10 34 -10 36 erg s _1 (Mereghetti 2000; Thompson 
2001). AXP sources likely have somewhat larger active ages 
and some of them have softer X-ray spectra compared with 
SGRs. One of the primary difference between them will be 
that AXPs so far have shown only quiescent X-ray emission 
with no bright active bursts such as giant flares. Active ages 
of SGRs and AXPs are consistent with the observed evidence 
that these compact objects often give their location close to 
the edge of shell-type supernova remnants. Neither SGRs 
nor AXPs show the presence of conspicuous counterparts at 
other wavelengths (Mereghetti & Stella 1995; Mereghetti et 
al. 2002). Origin of these enourmous energetics is likely to 
come from their strong magnetic fields B = 10 14 -10 15 G es- 
timated by their spin periods and spin down rates. Only the 
magnetar model has been able to account for the enigmatic 
properties of a rare class of SGRs or AXPs. 

Shear and Alfvenic waves play an essential role on the 
energy transfer to the exterior in the burst-like phenomena 
observed in SGRs or giant flares. The Alfven wave propa- 
gation in such a strongly magnetized star should be clari- 
fied theoretically. In our previous work, the propagation and 
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transmission of torsional Alfven waves have been studied, 
focused only on the fundamental mode of azimuthal wave 
number m = 1 as a first step. In that work, the exterior 
of the star is examined qualitatively only, assuming two ex- 
treme cases: (i) corotating together with the star and (ii) 
static state independent of the stellar rotation. The former 
will come true when the plasma gas is trapped by closed 
magnetic field lines, and the latter will be realized without 
any other force. In any case, it might be inadequate at least 
in the following two points to place such strong constraints 
on the wave mode and on ambient circumstances of the star 
in the previous paper. First, under the more realistic situ- 
ations various modes with azimuthal wave number m 3> 1 
are probably triggered by star quake. Second, it is natural 
to assume that a scale height of the plasma gas is neither 
L — nor L — > oo, but has a finite value of L. Huang et 
al. (1998) pointed out that a relativistic fireball like those 
in classical GRBs may exist in SGRs. Recently, Thompson 
and Duncan (2001) have also suggested that after the initial 
hard spike emission, some lumps of hot plasma gas involv- 
ing electron-positron pairs and high energy photons, that 
is, fireball, would be created on closed magnetic field lines. 
In fact, the fast decline and complete evaporation on X-ray 
light curve observed in the August 27 burst provides a clear 
evidence of the trapped fireball. Motivated by this, we thus 
extend our study to include some plasma gas spreading over 
the stellar surface. 

The main aim of this paper is to study the physical be- 
havior of the torsional Alfven waves with azimuthal wave 
number m in the presence of corotating plasma confined 
within a certain finite distance L and then to investigate the 
effects of these quantities m and L on the wave propagation 
and transmission. In section 2, our model is self-consistently 
constructed and the relevant basic equations are formulated. 
These equations are almost the same as those derived in 
Kojima and Okita (2004), but are summarized here for the 
paper to be self-contained. In section 3, we shall derive the 
dispersion relation in a WKB way to show that the rotating 
plasma atmosphere plays a crucial role as a resonant cavity 
of the wave when a certain condition holds. In section 4, 
transmission rates of the Alfven waves are numerically cal- 
culated. In section 5, their electromagnetic field structure is 
also discussed with the numerical results. In section 6, we 
give a brief summary of our findings and their implications 
for the torsional Alfven waves on rotating magnetars. 



2 ELECTROMAGNETICS ON ROTATING 
MAGNETARS 

2.1 Model 

Both the magnetic field and the rotation of star lead to a 
complicated geometrical configuration. In this paper we as- 
sume some simplified conditions to understand the physical 
processes of the wave propagation. We here consider three- 
layered cylindrical model with radius vj pc - It is composed of 
the neutron star crust between z = — q and z — (region 
(1)), corotating plasma above the stellar surface between 
z — and z = L (region (2)) and static pure vacuum at 
z > L (region (3)) as shown in Fig.l. The plasma gas filled 
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Figure 1. Three-layered cylindrical model. Region (1) corre- 
sponds to the rotating neutron star crust with a constant angular 
velocity Q, region (2) corresponds to the fireball as a resonant cav- 
ity filled with plasma corotating with S7, and region (3) denotes 
a static pure vacuum. 



in the atmosphere corotates with the crust at the same an- 
gular velocity Q. Local magnetic fields permeate uniformly 
in each layer and point along the z direction, B = B e z , 
which represents open magnetic fields extending to infinity. 

Alfven waves excited at the bottom of the crust q « 10 5 
cm owing to some mechanisms, travel upward along the local 
magnetic field lines. They are partially reflected and trans- 
mitted at the boundaries z — and z = L, where the physi- 
cal property of Alfven waves significantly changes because of 
the effect of the background rotation, as shown below. This 
model, by setting L — > oo, reduces to the case in which the 
plasma gas extends infinitely to the exterior of the star, while 
setting L — > reduces to the simple case in which the exte- 
rior is filled with pure vacuum only. Up to the present time, 
we have no observational information of the plasma size L. In 
this work, L is regarded as a free parameter in order to inves- 
tigate its effect on the wave propagation and transmission. 
We shall now give some comments on the validity of this 
model by comparing the physical sizes vj pc , q and L with the 
stellar radius R. In this model we have explicitly assumed 
that the local magnetic field has a ^-component only. Our 
model can be applied to the polar cap region, whose cylindri- 

-1/2 



cal radius zo pc is given by zu pc 



Rsme pc « 10 4 (T/ls) 



cm <JC R — 10 6 cm. Curvature of the stellar surface and 
the magnetic field lines may be neglected within the polar 
cap region. In a similar way, it may be valid to assume that 
the local magnetic field lines are uniform if the thickness of 
each layer is smaller than the star radius, L,q <C R- We can 
also extend our model to the extreme case L ~ R. Even in 
this case, the validity of this model nevertheless holds good 
near the z-axis. In the remainder of this paper we ristrict 
our explicit calculations only to the axially symmetric small 
region within the polar cap radius, unless otherwise stated. 
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2.2 Linear Perturbation 

In this section, we consider the propagation of torsional 
shear-Alfven waves with various azimuthal modes on rotat- 
ing magnetars. Such waves are probably excited by the tur- 
bulent motion of the starquake in the deep crust, but above 
the neutron drip (z > z n d ~ — 10 s cm). Many proposals 
for starquake model have been put forward (e.g., Pacini and 
Ruderman 1974), but all of them are generally argued only 



for weakly magnetized neutron stars with 10 



10 1 



Therefore, their treatments are inadequate for magnetars, 
as they are. However, analogous mechanism may also occur 
in magnetars. In this paper, we do not discuss the trigger- 
ing mechanism of the starquake itself, but focus only on the 
process by which electromagnetic shear waves, once gener- 
ated, are propagated and transmitted from the deep crust, 
through a magnetized plasma, toward the vacuum exterior. 

We assume that the horizontal Lagrange displacement 
£ = (£ ro ,£0,O) is suddenly shaken in the deep interior at 
depth q ~ 10 5 cm, nevertheless the matter remains immobile 
in the vertical direction due to strong gravity of the neutron 
star, g ~ 10 14 cm s -2 . Therefore the following transverse 
wave condition may be easily satisfied. 

1 



1 d 



= 0. 



(1) 



As one of the simplest solutions, wave ansatz under this 
condition can be written formally as 

m — 1 



(2) 



i{uj,4>, z,t) = (e ro ± ie^) ( ^— ) f±m,w (2) 
(m = l,2,3---)- 



Here notation '±' denotes the difference of helicity states. 
Displacements and other perturbed quantities are always 
assumed to have a harmonic time dependency e -1 "*. 

The crust and the surrounding plasma can be regarded 
as a perfect conductor because ohmic dissipation timescales 
are very long compared with other timescales of interest. 
The frozen condition can read 



E+-v x B =0. 

c 



(3) 



Background is assumed to rotate uniformly with v = tufle^, 
which means that the local electric field arround a highly 
conducting spherical star has only a radial component 
E = (— wQ,B /c)e m . The electric field now generates the 
Goldreich- Julian charge density 



p c = V • E/4tt = ~nB /2nc. 



(4) 



Transverse displacements of the crustal matter produce 
electromagnetic field perturbations propagating along the 
local magnetic field B = B e z . A self-consistent set of 
Maxwell's equations for perturbed electromagnetic fields <5E 
and SB arc 



V • <5E = 4tt8 p , 
V ■ SB = 0, 
1 



V x 5E : 



dtSB, 



V x <5B = —8] + -d t 8E, 

c c 



(5) 
(6) 

(7) 
(8) 



where S denotes Eulerian perturbation. By taking the per- 
turbation of equation © and then combining with equation 



Q, <5E and <5B are expressed as 
1 



SE - 



(Sv x B + v x SB) , 



SB = (B V)£. 



(9) 
(10) 



By using the relation between the displacement and the ve- 
locity perturbation Sv = <9t£ + (v • V) £ — (£ ■ V) v together 
with equation 0, one obtains 



SE=^ 



=F (u! =p mO,) f± (e ro ± ies) + e z 

dz 



SB = Bo I — z— (e ro ± ies) e v ^ 

, vjpc J dz 



(11) 



(12) 



Equations 1111 and (1121 1 imply that if stellar rotation can be 
completely ignored, then SE, SB and the propagation vec- 
tor k = ke z form a mutually orthogonal set of vectors, say, 
TEM modes. However, in the presence of rotation, such an 
orthogonality breaks down and thus longitudinal modes of 
the perturbed electric fields and currents are excited (TM 
modes). Other perturbed quantities related to the displace- 
ment are calculated by using above the results (II IP and (1121 



8j = — V x SB 



-d t SE 



_.CB 



4tt 



4-k 

m—1 



dz 2 c 2 

— i(ujt^fm<f>) 



£± > (e ro ± ie^) =p 



loQvj d£± 
c 2 dz 



8 Pc = —V ■ SE 

47T 

B o 0,lu ( vj 



4irc 



m—l ,2 C 

p. 

dz 2 



(13) 



(14) 



The linearized equation of motion for uniformly rotating 
background from a static equilibrium state governing the 
torsional Alfven waves is now given by 



p [dtSvi + (v • V) Sin + (Sv ■ V) Vi ] 
SS, + SFi + ( S Sp) i - (VSp). (i = vo, 



(15) 



Profile of the mass density p in the crust is obtained by solv- 
ing the equation of state for degenerate electrons as follows 
(Blaes et al. 1989) 



(p e m u ] 



5/2 



3n 2 h 3 
8.0 x 10 



g p c m u 2 , „ |, 
z +2gm e \z\ 



+ 2.5 x 10" 



3/2 



3/2 



g cm 



(16) 

where m u is an atomic mass unit, p e is the mean molecular 
weight per electron and others have usual physical meanings. 

The first term of RHS in equation 115H . 8 Si, is a per- 
turbed elastic stress tensor associated with distorted matter 
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in the crust and can be written in terms of £ 



SSi = V, 



dc_ + oe_ 



= m,0,z), (17) 
where x is a bulk modulus, is a shear modulus given by 



^ = 0.295Z 2 e 2 n* /3 



4.8 x 1(T 



10 11 g cm- 



4/3 



erg cm 



(18) 



with the ion number density m = p/Zp c m u and <5ij Kro- 
necker's delta. In the following calculations, Z = 32 will 
be adopted as a typical value in the crust. If one considers 
the complete transverse oscillation mode, the first term of 
equation 1171 vanishes owing to V • £ = 0. 

The second term on RHS of equation (1151 can be re- 
solved into perturbations of Lorentz and Coulomb forces 



5F = 5p c E+p c SE+- (Sj x B + j x SB) 

c 

— ~ (<5j-pe5v) X B. 



(19) 



Here we used equation © and unperturbed current j =p c v 
induced by the background rotation. We further dropped 
the term 8p B E oc (fizu/c) 2 , which is small near the z-axis or 
within the actual stellar radius R. Eliminating 5j with the 
help of equation JHJ, the electromagnetic force reduces to 



*F~_CB.V)'f. 



47TC 2 



(20) 



The first term in equation (12(H implies the tension of the 
perturbed magnetic field, while the second one shows the 
magnetic pressure generated by the distorted matter. 

The last two terms of equation l|15^ represent gravi- 
tational force (gSp) i and pressure gradient (X75p) i , respec- 
tively. However, unless one takes the p- or f-mode such as 
compressional waves and/or sonic waves into consideration, 
one can ignore them for simplicity in this model. 



2.3 Wave Equation 

We shall now derive the wave equation in region (1) shown 
in Fig.l. Substituting equations 11711 and 12011 into equation 
of motion 115H . one can obtain the following differential wave 



equation in terms of the displacement £ 



(i) 



d 2 ^ 1 dp dZ 



Ci) 



dz 2 



■ + 



p dz dz p 



+ ^a ± [a± ± {m(l - h) + 2} ft] $ = 0, 



(21) 

where p and p denote the effective shear modulus and the 
effective mass density defined as 



p = p + 



Bl 

47TC 2 ' 



(22) 
(23) 



The ratio of these quantities gives the shear-Alfven wave 
speed in the crust v = ^/p/p. In the above expression, the 
frequency cr± = uo =p mQ. measured in corotating frame for 
each helicity state has been introduced. We can limit this 
frequency to the positive regime a± > for the symmetry 



£±m,-w = £:pm,w Dimensionless function h in equation 12H 
is formally defined as 



h ■ 



4npc 



Anpc 2 + Bl 



(24) 



which has a great influence on the dispersion relation of 
the wave especially with large m not only in the inner sur- 
face, but also in the rotating plasma cloud, as discussed in 
the following section. Note that if one considers the static 
background (fl = 0), wave equation 12 U coincides with the 
one already derived by Blaes el al. (1989). We now look for 
a WKB solution. In the deep interior the solution of wave 
equation (12 1 1 can be well asymptotically given by 



^asymp 



\zf {A± exp [-i(ip±(z) +ut)] 



+B±exp[i(1>±(z)-vt)]}, (25) 
with /3 — —7/4. Here tp±(z) denotes eikonals defined as 



ip±(z) 



dz * ^/a ± [*± ± {mjl - h) + 2} H] 



for each mode. The first and second terms in equation l|25fl 
represent the upward-propagating Alfven wave with a com- 
plex incident amplitude A± and a downward-propagating 
wave with a complex reflection amplitude B± bounced at 
the stellar surface, respectively. 

We now turn to the wave behavior in region (2) (0 < 
z < L). Mass density in this region is so small that one can 
formally take the limit h —> 0. Thus equation 12 U reduces 
to 



dz' 2 



+ ^<J± [<r± ± (m + 2)Q] Ci 2) = 



(27) 



The solution of this equation can be analytically written as 

££ 2) = C± exp [i(k (2) z ~ ut)] + D± exp [~i{k (2) z + ut)] , 

(28) 

12) 

where the wave number k ± with each mode in the plasma 
is defined as 



-V°± k± ± (m + 2)n]. 



(29) 



In region (3) with pure vacuum (z > L), the wave equa- 
tion becomes 



— 1 £ ( ' = 0. 

dz 2 c 2 ? 



(30) 



Owing to the absence of rotating matter, two helical states 
satisfy the same equation. We here dropped the notation 
'±'. The solution is easily given by 

£ (3) = Eexp \i(k {3) z - ut)] ■ (31) 

The wave number thus takes an ordinal form as 

fc (3) = - (32) 



c 



In the above expressions, C± , D± and E are complex inci- 
dence, reflection and transmission amplitudes in each region, 
respectively. These wave amplitudes and the wave numbers 
determine the transmission rate of the wave based on some 
boundary conditions. Mathematical treatment will be given 
in section 2.4. 
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2.4 Boundary Condition 

The physical property of Alfven waves changes at the bot- 
tom and at the top of the rotating plasma layer, depending 
on the azimuthal wave number and the angular velocity of 
the background. We now require boundary conditions in the 
usual way in order to connect each wave solution continu- 
ously. From equations 1281 and I31H . the continuity at the 
upper surface of the plasma, z = L, gives 



u(2) 



c± 



E- 



D± — E 



2k 

fci 2) 



(2) 



■exp 



.(2) 



u(3) 



L 



,(3) 



2k 



(2) 



■ exp 



[i (4 2 



+ k 



(3) 



L 



(33) 



(34) 



Substituting equations l|33|l and l|34fl into equation l|28fl and 
then taking the derivative at the stellar surface, one obtains 



— In 

dz 



t(2) 



(35) 



where 7 is a modification quantity due to background rota- 
tion of the surrounding plasma denned as 



7 = 



>c+ exp [- 



exp [ix+L] 



x+ exp [—i>c-L\ + x- exp [ijt+L] ' 



with 



= *4 2) + k 



= k 



(2) 



(3) 
k^. 



(36) 
(37) 



Detailed treatment and physical meaning of the 7 are given 
in the Appendix. 

In order to solve wave equation 12 U . we consider the 
complex linear combination C — Ci + 2C2 with specific so- 
lutions fi and £2 to be satisfied with equation 1351 . At the 
ends, one can obtain the explicit boundary condition at the 
stellar surface z = 0, 



d_ 

dz 

= -fci 2) {Ci(0p 7 + C2(0)5R 7 } 
+ikf {£(0)3*7- <2(0)$>y} • 



(38) 



Logarithmic derivative of asymptotic solution has to be 
equal to that of numerical solution in the deep interior 
2 <C 0, so that we request 



d_ 

dz 



In (£^ m P) 



= E ta <o 



This yields 
B± = _ 



-1 . .dt/}± 
dz 



c + 



dC 

dz 



-1 , { di/j± 
dz 



^ dz 



exp [— 2iip±] 



(39) 



(40) 



Finally, the reflection R± and transmission coefficients T± 
of the waves with each helicity propagating from the crust 
through the plasma toward the vacuum exterior are respec- 
tively expressed as 



R± = 



\B±\ 



(41) 



k& \Ef 
fc« \A ± 



R±. 



(42) 



3 PROPAGATION 
3.1 Dispersion Relation 

In this section, behavior of the torsional Alfven waves is dis- 
cussed based on the dispersion relations. Eikonal equation 
1261 depends on the depth through the functions h and v. 
For strong magnetic fields B ^ 10 1 G considered in this pa- 
per, it can be shown that tjj± in equation 1261 varies weakly 
with the depth except for the inner region close to the star 
surface. The local dispersion relation can be approximately 
written as 



fc± W - \/o± [a± ± {m(l - h) + 2} ty . 



(43) 



This relation can also be obtained by taking a short wave- 
length limit, namely, high frequency limit. The dimension- 
less function h involved in equation 1)43^ varies from h = 1 
in the deep interior to ft = in the surface and exterior. 
This function represents the ratio between the rest mass en- 
ergy density and the effective magnetic energy density. For 
h — 1, the rest mass energy dominates over the magnetic 
energy. This case corresponds to the classical limit. On the 
other hand, strong magnetic energy B a ^ 10 14 G easily over- 
whelms the rest mass energy of the electrons or ions in the 
low density region, where h ~ 0. In this case, relativistic 
displacement current should be included in the analysis. 

In order to investigate the propagation property of the 
wave, it is usefull to define the phase velocity v p in a straight- 
forward manner as 



u p/c 2 = '/ (c/c±) 2 



(44) 



which is also related to the refraction index N± = c/v p 
for each mode. Note that Vp is not necessarily positive. It 
is well-known in plasma physics (e.g., Wolfgang and Rudolf 
1996) that waves are in general reflected at the cut-off points 
where N± — > 0, and are absorbed at the resonant absorp- 
tion points where N± — > 00. Since the frequency a± in the 
rotating frame is here confined within the positive regime 
a± > 0, fc+ is always real, that is, neither cut-off nor res- 
onant absorption points appear in the positive helicity. On 
the other hand, there are both absorption and cut-off points 
for the negative mode. In the remainder of this section, our 
attention will be thus paid only to this mode for physical 
interests. The dispersion relations with negative mode are 
schematically shown for classical limit in Fig. 2 and for rel- 
ativistic limit in Fig. 3, respectively. In both diagrams the 
region fp/c 2 < stands for the non-propagation, viz. evanes- 
cent zone. Vertical dotted line denotes the cut-off frequency 
of the wave. 

(i) Classical Limit 

As seen in Fig.2, evanescent zone appears only in the 
low frequency and narrow band width < cr_/f2 < 2. In 
this classical treatment, the cut-off of the wave appears at 
= 2, which can be physically interpreted as the fact 
that Coriolis force interrupts the wave propagation in the 
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Figure 2. Dispersion relation for negative helicity waves with 
m ^ 3 in a classical limit, h — > 1. The region below horizontal 
axis corresponds to the evanescent zone. Vertical dotted line de- 
notes the cut-off frequency cr_/n = 2. The resonant absorption 
frequency is given by <j—/il = 2. 
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Figure 3. Same as Fig. 2 in a relativistic limit, h — ► 0. The cut- 
off frequency is cr_/Q = m + 2. The resonant absorption occurs 
within an evanescent regime at cr_/f2 = m. relation for negative 
helicity waves with m ^ 3 in a classical limit, h — > 1. The region 
below horizontal axis corresponds to the evanescent zone. Vertical 
dotted line denotes the cut-off frequency <r_/f! = 2. The resonant 
absorption frequency is given by cr_ /Q = 2. 



rotating frame. In a high frequency region a-/Q, > 2, waves 
are almost capable of propagating except for <r_/f2 — m, 
at which waves are strongly absorbed by resonance with a 
rotating background. This is because refractivity diverges in 
this particular frequency. One notices that such a resonant 
absorption frequency <j_ = mQ, in our model is formally 
analogous to the electron-cyclotron resonant frequency in 
the standard plasma physics (Wolfgang and Rudolf 1996). 

(ii) Relativistic Limit 

It is important to understand the new effect that ap- 
pears in the relativistic case. Figure 3 demonstrates that 
evanescent zone prevails in the frequency regime < 
<r_/n < m + 2, whose band width is therefore broadened by 
large m modes. Cut-off frequency is given by cr_/Q = m + 2. 
It is only a high frequency wave with <r_/fi > m + 2 that 
can always propagate in a WKB sense. By contrast with 
the classical limit, resonant absorption a-/fl = m does not 
appear in this propagation regime, but in the evanescent 
regime. 




Figure 4. Spatial structure of the potential barrier for the nega- 
tive mode as a function of distance. Horizontal dotted line denotes 
rotation frequency cr_ of the wave propagating above the poten- 
tial. 



3.2 Structure of the Potential Barrier 

As is obvious from Fig. 3, evanescency can not be neglected 
in large m mode. In this subsection, we turn to the inves- 
tigation of the spatial structure in the evanescent zone for 
large m. We now introduce a new dependent variable ^_ 
defined by ^_ = ^ 1//2 £_ for the amplitude of the negative 
mode with z = y/a-z. Exploiting some arithmetic algebra, 
wave equation 112 It is now rewritten into a Sturm-Liouville 
type differential equation 

\ rr" a i / ■ 

Vrr. 



dz 2 



+ 



4/+ 



( *hL 

\ dz 



where 



2ju dz 2 



E : 



9- 



= 0, 
(45) 

(46) 



Vrr, 



= {m(l-h)+2}Sl. (47) 

The first two terms concerning Jl in square brackets in 
equation I45H become very small for strong magnetic fields 
B > 10 14 G, since Jl can be almost regarded as a constant. 
One can compare equation 1451 with the one-dimensional 
Schrodinger equation for box-type potential in a stationary 
state. Physical quantities E and V m respectively correspond 
to wave energy and potential in a formal sense. Strictly, po- 
tential Vm depends on the wave energy E through the fre- 
quency <t_. However, cr_ is here treated as a constant value 
irrespective of the position, so that the following discussions 
are valid without loss of generality. 

Schematic profile of the effective potential V m and the 
wave energy E are given in Fig. 4 as a function of distance 
z. We draw the curve in region (1) somewhat exaggeratedly. 
The potential varies with the position through the function 
h and the background rotation Q,. The potential height is de- 
termined by the coupled quantity (m + 2) Q and its width 
is given by the size L of the corotating zone. This result 
means that the large wave number lifts up the potential bar- 
rier, only if the background stellar medium rotates. In other 
words, if surrounding matter is not dragged by the star, say, 
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being in static state Q — 0, then the potential never rises, 
even though some high azimuthal waves exist. 

Whether or not the wave can propagate and transmit 
out depends on the wave energy, that is, wave frequency 
<r_. Our argument deserves to be specially emphasized in 
two explicit regimes; the evanescent frequency mode 2Q < 
a- < (m + 2)Q and the propagation mode er_ > (m + 2)Q. 

3.2.1 Evanescent Mode: 2fi < a_ < (m + 2)Q 

As is apparent in Fig. 4, such low frequency waves excited in 
the interior cannot help striking the potential barrier. Since 
refractivity is zero on the critical curve cr_/f2 = m (1 — h) + 
2, outgoing waves with m (1 — h) + 2 < cr_/f2 < m + 2 are 
generally reflected, when they reach the potential wall. 

3.2.2 Propagation Mode: cr_ > (m + 2)Q. 

In this case perturbation of the electromagnetic fields can 
propagate as a wave. The present context in our model is 
almost concerned with the WKB frequency range. The prop- 
agation sometimes exhibits a remarkable property, if certain 
conditions are satisfied. As long as the wave frequency is 
much greater than the potential strength cr_ 3> (m + 2)Q, 
the potential itself does not have much influence on the wave 
behavior. However, if the frequency becomes commensurable 
to the potential height <r_ > (m + 2)f2, the existence of a 
potential barrier cannot be ignored. Especially if the wave- 
length is comparable to the potential width, i.e., k-L ~ I, 
then the waves will interfere with the potential barrier and 
their behavior will be strongly altered. This inherent prop- 
erty is expected to be more evident near the threshold fre- 
quency cr_ ~ (m + 2)0,. Motivated by this general con- 
sideration, we explicitly calculate the transmission rate in 
the next section. The numerical parameters are adopted 
to satisfy the above condition, L ~ l/fe_ ~ 10 6 cm and 
m = Ld max /fi ~ 10 6 . Using these parameters, we will exam- 
ine how and to what extent the propagation and transmis- 
sion are affected by the potential barrier. 



4 TRANSMISSION 

4.1 Numerical Calculation 

We numerically calculate the transmission rates (1421 of the 
wave propagating from the deep crust, through the coro- 
tating plasma envelope, toward the vacuum exterior sub- 
ject to matching conditions at each boundaries. As already 
mentioned above, only the negative helicity wave is intrigu- 
ing for physical interests. Evanescent modes in a low fre- 
quency regime < er_ < {m(l — h) + 2} Q are excluded 
from this calculation. This constraint on wave frequency 
thereby guarantees that all waves are capable of propagating 
in a WKB sense. One can roughly estimate a typical wave 
frequency ui measured in the inertial frame by approximat- 
ing as ui ~ v/q. Appropriate frequency thus lies in a finite 
range I0 3 < w < I0 6 s" 1 . 

The objective of this section is to explore the effects of 
azimuthal wave number m and the size L of the corotating 
plasma medium on the wave transmission. We here consider 



two kinds of specific torsional waves whose azimuthal num- 
ber is; (i) the fundamental, m = 1 and (ii) much larger 
than one, m = aj m a X /0 = I0 6 . In each case, we further con- 
sider two apparently different circumstances in the exterior 
L = R ~ A_ ~ 10 6 cm and in the absence of plasma L = 
cm. 

Figures 5 and 6 respectively demonstrate the transmis- 
sion coefficients for the negative helicity waves with two 
extreme cases (i) m = 1 and (ii) m — 10 6 as a function 
of the wave frequency measured in the inertial frame when 
B = I0 15 G, typical for magnetars. In both figures, solid 
and dotted lines denote the results of L — 10 6 and L — 
cm, respectively. The surrounding plasma is here assumed 
to corotate with the same angular velocity as that of the 
star, f2 = 1 s _1 . In each case, we got the following results. 

(i) The Fundamental Wave Number: m = 1 

As seen in Fig. 5, the transmission curve in the case of 
L = I0 6 cm slightly shows a wiggling behavior. This curve 
intersects that of L — at u ~ 10 4 s _1 . However, the 
difference between both cases becomes undistinguishable in 
the high frequency regime u > 10 5 s _1 , since the wave fre- 
quency is much higher than the critical one (m + 2) fl — 3 
s _1 , which appeared as the potential barrier in section 3. 
As the wave frequency gets higher, the transmission rate 
approaches unity asymptotically. Since the overall property 
does not depend on L, the potential barrier due to rotating 
plasma does not make any significant influence on the wave 
transmission for small azimuthal number m ~ 1. 

(ii) Highly Azimuthal Wave Number: m = I0 6 

In this high m mode, we can find some suprising results. 
Figure 6 shows that the transmission rate of L = 10 6 cm 
is drastically enhanced at some selected frequencies. More 
importantly, such enhancements occur periodically at spe- 
cific frequencies. At the top of the first wing (n = 1) u)% = 
2.2 x 10 s s _1 , the rate reaches the maximum T max =0.7I, 
which corresponds to approximately 70 times larger than 
that of the first bottom ui D = 7.0 x I0 3 s _1 . Comparing with 
L = case at the same frequency uj t , this maximum T max 
amounts to 230 times larger. In this way, plasma effect is 
important at each top frequency, but unimportant at any 
bottom frequency. As shown in Fig. 6, the transmission co- 
efficient through the plasma layer is exactly equal to that 
of L — case at arbitrary wing bottom. In general rotating 
plasma has an effect in helping the escape of the waves ex- 
cept for the bottom frequencies. This point is different from 
that of case (i). 

In order to understand the transmission enhancement, 
we have also numerically computed the integral J Tdui in 
some frequency bands for both L = I0 6 and L — 
cm. We obtain some 13 times enhancement in the nar- 
row band 2Q < uj < uj a and 1.4 times in the broad band 
20 < u < w m ax compared to the results of L = 0. Such 
enhancement and periodic variation slowly decrease with 
increasing wave frequency. This property almost vanishes 
when the frequency becomes comparable to the critical fre- 
quency (m + 2) Q, ~ 10 6 s" 1 . 

Periodic enhancements in our numerical calculations 
can be interpreted as a consequence of wave interference 
within the rotating plasma cavity. In the following subsec- 
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this homogeneous model, the transmission coefficient T can 
be analytically calculated by taking boundary conditions at 
the discontinuities z — and z — L, 



T 



+ 



4fc' 3 > 
f fcW 

\ fe( 2 ) 



l+m] cos 2 (fc< 2 >L) 



fe(i) 



2 (fc (2) L) 



(48) 



Although the physical situation is different from our 
present work, similar expressions could be found also in our 
model. This formula 1481 implies that the transmission rate 
has a periodic structure depending on the wave number 
and the characteristic scale L of the intermediate layer unless 
k^L ~ 0. Extrema of the transmission rate work out to be 



Figure 5. Transmission coefficients of the negative helicity wave 
with the fundamental mode m = las a function of frequency when 
-B = 10 15 G and f2 = 1 s —1 . Solid and dotted lines denote the 
results of L = 10 6 and L = 0cm, respectively. 




Figure 6. Same as in Figure 5, but for high azimuthal number 
m = 10 6 . In this high m mode transmission rate is highly en- 
hanced at some selected frequencies due to the plasma envelope. 



tion, we investigate our results more quantitatively in a sim- 
plified analytical method. 



4.2 Analytical Calculation 

Periodic enhancements on the wave transmission in our 
model can be satisfactorily rationalized by comparing with 
the wave propagation in homogeneous multi-media. Non- 
relativistic Alfven resonance itself in the magnetosphere 
around the earth and the sun have been widely recognized 
and discussed theoretically and observationally (for recent 
reviews see Leonovich and Mazur 1997; Waters 2000; Li and 
Wang 2001). Hollweg (1983) has theoretically studied WKB 
wave propagation in three homogeneous layers labeled by 
(1), (2) and (3) separated at z = and z = L as shown in 
Fig. 7. In his work Alfven waves are simply assumed to have 
constant wave numbers k^\ fc' 2 ' and k^ in each region. In 



?max — 



4fc< 1 >fc( 3 > 

(feCl) +fc(3)) 2 



at fc (2) L = n-K, 



4fc( 1 'fc( 3 V(fc( 2 >) 2 

(l + ft(l)fc(3)/(fc(2))2) 



at k {2) L = (2ra- l)7r/2, 
(49) 

with n — 1, 2, • ■ • . Providing that the wave number in each 
region satisfies the inequalities 

fc (1) »fc (2) »fc (3) , 

k (2) L > 1, 



(50) 
(51) 



the coefficient of cos 2 (fc ( - 2 - ) L) in equation 14811 becomes domi- 
nant compared to that of sin 2 (fc ( - 2 'L). In this limit, equation 
can be well approximated by 

4fc< 1 >fc< 3 > 



(fc (2) L). 



(52) 



(Jfc(i) +fc< 3 )) 2 

One notices that if k^ ~ k^ ~ fc' 3 - 1 , the transmission 
shows neither periodic variations nor enhancements. Conse- 
quently, inequalities J5UI and (I5H give a set of resonant con- 
ditions of the wave. Sterling and Hollweg (1984) have subse- 
quently considered a three layer model for solar flare, which 
is composed of a chromosphere, a spicule and a corona. In 
that work they have confidently suggested the possibility of 
Alfvenic resonance on solar spicules and have shown a new 
aspect of spicules which may account for occasionally twist- 
ing motions of magnetic field lines, even when above res- 
onant conditions approximately hold. Conditions (46) and 
(47) may hold true also in our model except close to the thin 
regime beneath the stellar surface, whenever large twisted 
waves with m^> 1 propagate in the rotating background. 

We can now quantify the particular frequencies at which 
the wave resonance occur. By solving the quadratic equation 



(<T- 



(m + 2) ricr_ 

(2), 



( 2 h 



riodic condition k_L= (2n 



= 0, together with the pe- 
1) 7T, the eigen frequencies of 



the negative mode are obtained as 



(2 



-fi + 



m + 2. 



n + 



(2n - 1) C7T 
2L 

' (53) 

with m,n — 1, 2, 3, • • • . This formula shows that the reso- 
nant frequencies depend critically on the azimuthal number 
m and the plasma cavity length L. They gently decrease 
with an increase in L or m. Physically, this means that it 
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(1) 


C2) 







L 



Figure 7. Three homogeneous model separated by two discon- 
tinuities at z = and z = L. If periodic boundary condition 
= n7r (n = 1,2, •■•) holds, intermediate layer acts as a 
resonant cavity and transmission is highly enhanced. 



takes longer time for the wave to go back and forth between 
the stellar surface and the top of the plasma layer and then 
this wave interferes with another one propagating from the 
crust into the cavity. In extremal cases, one finds 



lim u) m ,n = hm u) m ,n = 2£1, 

L — >oo m — >oo 



(54) 



which coincides with the cut-off frequency in the classical 
limit. At the same time, transmission peaks become blended 
with neighboring resonance, 



lim Au) m — lim Acu m — 0, 

L — >oo m — >oo 



(55) 



-i- 



where Auj m = L> m ,n — Wm,m- 

For Q — 1 s _1 and L — 10 6 cm, from equation 15311 

one can calculate some representative resonant frequencies; 

8.8 x 10 3 s~\ 

and the third 

Our numerical work also gives sim- 
^S „— 1 , . i n w i n4 ., — 1 



the fundamental mode frequency u 10 6 1 - 
the second mode lo 10 6 >2 = 7.4 x 10 4 s" 
w 10 6 i3 = 1.9 x 10 5 

ilar results; tu 10 e tl = 2.2 x 10 J s" 1 , uj 10 6 t2 = 1-9 x 10* s 
and uj 10 6 3 = 5.2 x 10 4 s _1 , respectively. Our results exhibit 
slightly positive deviations from the analytical ones. Recall 
that the wave number in our magnetar model cannot 
be regarded as a constant, since the mass density drastically 
changes in the vicinity of the stellar surface. This discrep- 
ancy would be probably attributed to the inhomogeneity in 
the crust. If a sharp boundary is formed at the stellar sur- 
face, such a difference may probably become small. From 
a comparison with the analytic model, we concluded that 
the transmission enhancements obtained in our model are 
thought to be a result of the Alfven resonance on the rotat- 
ing plasma cavity due to the potential barrier. 



5 ELECTROMAGNETIC FIELD STRUCTURE 

In the preceding section we have elucidated that the trans- 
mission rates of wave are highly enhanced at some selected 
frequencies owing to the resonance effect in the rotating 
plasma cavity. Resonance in the plasma portion has another 
spectacular nature. In this section, it will be shown that 
resonance has a great impact not only on the transmission 
rates, but also on the electromagnetic field structure associ- 
ated with Alfven waves. Most of our applications will once 
again concern the WKB frequency regime of the negative 
helicity er_ > {m(l — h) + 2} fl. Hereafter we omit the sub- 
script ' — ' for simplicity. 

By substituting equation 1281 into 1111 and 111 we 
can explicitly have the expressions for electromagnetic field 
amplitudes in the plasma zone in terms of ft' 1 ', normalized 



by a local magnetic field B 



B- 



B 2 

ft, 



\cSE z 



Q, 2 uj 2 Bl 



(2) 



dz 



(1) 



1-41 



X [(feL 2) ) 2 sin 2 k {2) (z — L) + (fc (3) ) 2 cos 2 k {2) (z - L)J , 



(56) 



and 



\c8E m \ 2 \cSE^\ 2 2 I (2) | 



/m(ro) (<7 m ) 



I A 1 Tm X 



{k^ym 

[(fcL 2) ) 2 cos 2 k% (z - L) + (fc (3) ) 2 sin 2 k% (z - L)J 



(57) 



where f m (zu) is a radial profile given by fmfo) = 
(tu/tUpc) 2 '" 1-1 - 1 (m = 1, 2, ■ ■ • ) and the other notations have 
the same meaning as the previous ones. The longitudinal 
component of the perturbed magnetic fields is always zero 
SBz — 0, because the matter has been assumed to be ver- 
tically immobile in the present work. On the contrary, only 
if the background rotates, the electric fields have a longi- 
tudinal component whose structure is essentially identical 
with that of the magnetic fields. Thereby we only have to 
investigate the transverse components {zu, <f>) of the fields. 

For the sake of exploring the dependence of the az- 
imuthal number m on the field structure, we once more re- 
strict our discussions within two kinds of extreme cases: (i) 
the fundamental mode m — 1 and (ii) the large azimuthal 
wave number m — I (3> 1). In this limit, formulae (49) and 
(50) can be well approximated by 



\6B„ 



\SB^ 



\c8E z 



Bl Bl Q, 2 uj 2 Bl 

f m {zu)k^k^T m \A\ 2 form = l, 



fc(3) 



Tm L4 



2 sin 2 k {2) (z - L) 



for m = 1, 



(58) 



and 
\cSE a 



B 2 



\cSEtf, 
Bl 



(i) 



k m 

fm{zu)a 2 n -^rrT m \A\'' for 771=1, 



l fe(3) 



/m(ro)cr„ 



T m \A\ z cos 2 k (2) (z-L) 

for 777 = /. 



fc(3) 



(59) 



Here we have dropped some small terms coupled to (ft' 3 ') 2 , 
since inequality ft™' 3> k^ 3> ft' 3 ' holds for large m. Equa- 
tions 15811 and 159H have the implication that the fields are 
almost constant for small m modes, but are sinusoidally 
changed with z for large m. 

In principle, the absolute value of the incident wave 



10 T. Okita and Y. Kojima 





Figure 8. Spatial structure of perturbed magnetic field with 
m = 1 in the rotating plasma region normalized by that in the 
absence of plasma at some resonant frequencies in the first wing. 
The solid and the dotted lines denote the result at the top of 
wing ojt = 2.2 X 10 3 s _1 and at the bottom ui h = 7.0 X 10 3 s _1 , 
respectively. 



amplitude \A\ cannot be determined in our linearized the- 
ory. If we are allowed to assume that f m =i (zu) \A\^ =l ~ 
/ m -i (zu) |j4.|JL =1 , then the ratios of electromagnetic field 
amlitude with m = I to those with m = 1 are approximately 
given by 



\SB. 



2 

1 1 m — l 



\SB, 

and 
\SEi\ 



I 2 

lm=l 



\SE,\ 



, in 
i + — 

LO 



, in 
i + — 

LO 



! 2n 

LO 



1+2 



Tm=i . 2 , (2) 

sin «!/ 



(60) 



J""" 1 cos 2 kf£} (z — L) , 

im=l 



(61) 

with i = zu,(j>. These quantities are much larger than unity 
because of the extra factor (1 + m/ui) except for some 
special locations z = nir/km and z — (2n — 1)7T /2km 
(m,n — 1,2, •••), which correspond to the nodes of the 
standing Alfven wave. 

We will now examine the effect of the resonance in 
the rotating plasma cavity, whose thickness is given by 
L ~ R = 10 6 cm, on the spatial structure of the perturbed 
electromagnetic fields. Let us designate by a the amplitude 
of the perturbed magnetic field within the plasma normal- 
ized by that in the absence of the plasma 



a{z) 



\SB; 



2 

\L=R 



Tl=r c 2 
T L=0 uj 2 



',(2) , 



(62) 



The first term in square brackets is related to the devia- 
tion due to high torsional modes in the rotating background 
from the dispersion relation in a vacuum. The ratio a is 
plotted against the height from the stellar surface for (i) 
m = 1 in Fig. 8 and (ii) m = 10 6 in Fig. 9, respectively. Solid 
lines denote the numerical results at the top of the wing 



Figure 9. Same as in Fig. 8 for m = 10 . Field amplitudes are 
drastically and periodically enhanced for large m. This enlarge- 
ment is attributed to the resonance transmission and the conser- 
vation of energy flux as a multiplier effect. 



2.2 X 10 d 

n,3 



, while dotted lines denote the bottom 
Lo-t, = 7.0 x 10 J s _1 . The transmission rates at resonant fre- 
quencies obtained in section 3 are appropriately used in this 
calculation. 

(i) The Fundamental Wave Number: m = 1 

As shown in Fig. 8, in the case of m = 1, the wave am- 
plitude keeps one order of magnitude in the plasma. This re- 
sult agrees well with the fact that the dispersion relation for 
small m ~ 1 is almost equal to that of the electromagnetic 



fields in pure vacuum, (fc 



/c . The background ro- 



tation therefore does not affect this small m modes. 



(ii) Highly Azimuthal Wave Number: 



10 6 



Some significant differences can be found in this high 
mode. The perturbed magnetic field strength a for m — 10 6 
at the resonant frequency to t is drastically amplified up to 

10 5 atz< 10 5 cm. Then the field strength a is fallen 
down to a — Tl=r/Tl=o ~ 230 at the top of the plasma 
layer, in which the first node of the standing Alfven wave 
is formed. Even at the bottom frequency lo^ the strength a 
has about 10 2 at z ~ 5 x 10 5 cm. 

In the same way, the normalized amplitude of the per- 
turbed electric field is expressed by 



(3(z) EE 



\5Ei\ 



\SE, 

Tl=r 
Tl=o 



2 

L=0 



mn 



(k 



(2K 



(Mm ) — T 



'.(2) 



/ 1- \ 



(63) 



which is plotted in Figs. 10 and 11. When m = 10 , the field 
amplitude j3 of io t has almost the same magnitude 10 s as 
a near the stellar surface z < 10 4 cm. But near the top of 
the plasma z ~ 10 6 cm, in contrast to a, (3 has a maximum 
5 x 10 7 which corresponds to the loops of the standing wave. 
Approximately, f3 is much greater than a because of the 
extra term (1 + run/co) 2 ~ 10 6 in equation 1631 . 
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Figure 10. Perturbed electric fields with m = 1 in the rotating 
plasma at resonant frequencies in the first wing. Solid line and 
dotted line respectively correspond to u>t = 2.2 X 10 3 s — 1 and 
uj h = 7.0 x 10 3 s _1 . 
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Figure 11. Same as in Fig. 10 for m = 10 6 . 



Such a field amplification can also be explained by con- 
sidering a conservation of energy flux F ~ puu [£| 2 . The 
phase velocity v of the mode m > 1 in the plasma at- 
mosphere is much smaller than that in the vacuum, as is 
confirmed in equations 1431 and 1441 . Assuming the same 
value |£| at the stellar surface, which is irrelevant to L, the 
amplitude |£| (or. v~ 1//2 ) is enhanced in the rotating plasma 
region so as to compensate for the slowing down of the wave 
propagation. 



6 SUMMARY AND DISCUSSION 

We have studied the propagation and transmission of tor- 
sional Alfven waves along the rotation axis on magnetars 
by using three-layered cylindrical model. If the intermediate 
plasma layer rotates with angular velocity O, the propaga- 
tion property for largely twisted waves with m > to/Q is 
drastically modified. This middle zone can be regarded as a 
kind of a potential barrier, whose height is explicitly spec- 



ified by (m + 2)Q. Waves having an angular velocity of 2fi 
originates from Coriolis force in the classical limit, while an 
additional quantity mil comes from the relativistic effect. 
This therefore means that the potential is lifted up higher 
in the large m modes on the rotating background. It should 
be emphasized that this new finding is purely attributed 
to having incorporated the displacement current into the 
model. 

We have numerically computed the transmission rates 
of the torsional waves with large value m 3> 1 driven 
in the crust, through the corotating plasma with a finite 
size L ~ R, into the exterior. We found that the trans- 
missions are strongly enhanced for large m w oj/Q, at se- 
lected wave frequencies satisfied with a periodic condition 
kL ~ (2n — 1) 7r/2 (n = 1, 2, • • • ). Such a transmission en- 
hancement arises because the rotating plasma forms a res- 
onant cavity which traps the wave energy by virtue of the 
strong reflections occuring at the transition region at the 
bottom and at the top of the plasma layer. Efficient trans- 
missions due to resonance is also compatible with the fact 
that the reflection of waves at the stellar surface can almost 
completely be extinguished. Magnetospheric waveguide or 
resonance cavity can thus actually generate a set of coher- 
ent eigenmodes for high m. Resonance may be a signature 
of fundamental processes by which high order torsional os- 
cillations of waves, if there exists, can transfer much more 
energy out of the magnetars especially in the low frequency 
regime. 

It has been well-known in quantum mechanics that 
when the potential width is comparable to (or integer times) 
de Broglie wavelength of electrons L ~ nh/m a v, incident 
wave flux is bounded by finely tuning a phase relation with 
inversely propagating waves (e.g., Schiff 1949). Then res- 
onance is excited in the potential zone and the wave flux 
spends much of its time in the resonance cavity. Similar res- 
onance may occur also in our macroscopic model, if the char- 
acteristic scale of the evanescent zone is comparable to the 
wavelength of the Alfven wave. 

We have investigated the spatial structure of the per- 
turbed electromagnetic fields propagating in the rotating 
plasma, when standing Alfven waves are formed. Resonance 
in the plasma cavity for high m modes has another spec- 
tacular aspect. At some points magnetic fields of the Alfven 
waves are drastically amplified even up to a few 10 5 times as 
large as those in the absence of plasma or in the static exte- 
rior. This amplification can be straightforwardly explained 
by both the transmission enhancement due to resonance and 
the conservation of WKB energy flux. 

In this way, once the Alfven waves are resonantly ex- 
cited and transmissions are greatly enhanced, a substantial 
part of the wave energy will probably be transferred into 
the ambient charged particles such as electrons or positrons 
confined within the fire balls associated with the bursts on 
magnetars. The charged particles entwined around magnetic 
field lines will then be violently swayed. These disturbances 
may eventually produce high energy 7-ray or X-ray emis- 
sions from the magnetars. More importantly, if magnetic 
field lines at some altitudes are strongly distorted from 
equilibrium state by resonance and plasma fluid contracts 
into the surface, the field lines with antiparallel components 
will approach each other. Such geometry of field lines will 
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give rise to possible magnetic reconnection type events as 
Thompson and Duncan (2000) have suggested. 

The behavior of resonant MHD waves within the corona 
or spicules erupted from the solar surface have been studied 
so far. In fact, field line resonance has been found to occur 
on magnetic shells in the magnetosphere of the earth. But 
their physical treatment is inevitably restricted to the very 
weak magnetic fields. Resonant property found in our mag- 
netar model is very similar to that in solar astrophysics or 
planetary physics. We can expect that the same mechanism 
works also on the relativistic Alfven waves even in a different 
environment, that is, in an extreme circumstance accompa- 
nied by very strong magnetic fields such as the magnetars. 

As a final remark, we shall address on the angular ve- 
locity Q of the star and the azimuthal wave number m. In 
our model, the rotating background is throughout assumed 
to have a slow angular velocity f2 = 1 s _1 , which is actually 
observed in SGRs and AXPs. We concentrate our treatment 
on the torsional waves with a specific value m = 10 6 as a 
possible limit of high m. As shown in our work, the strik- 
ing behavior of the resonance necessarily requires high order 
torsional modes, m > w/(l s _1 ) ~ 10 3 -10 6 . It is not clear 
whether or not such modes with high wave number m real- 
istically exist on the magnetars. However, one should keep 
it in mind that the local dispersion relation of the Alfven 
waves are almost determined by the coupled quantity mQ 
except for the high frequency regime. This means that the 
physical behavior of the wave with high m = 10 3 -10 6 on the 
slowly rotating background Q = 1 s _1 is essentially equiva- 
lent to that of small m = 1-10 3 and fast rotation Q — 10 3 
s _1 . Neutron stars are often expected to have been a rapid 
rotator with fl ~ 10 3 s _1 at a star-born period. Our results 
would thus become much more effective especially on young 
magnetars. 

Appendix 

The physical effects of 7 can be clarified by taking 
some explicit limits. Specifically if one ignores the substan- 
tial thickness of the plasma layer L — *• 0, then the wave 
number of region (2) should be replaced with that of re- 
gion (3), *4 2) -> fc (3) . This yields that x+ -> 2fc (3) , h- -> 
and 7 — » 1. In this limit, the boundary condition reduces to 
d/dz[\n£±^] — ik^\ which clearly corresponds to the simple 
case that the exterior of the star is filled with pure vacuum 
or static plasma gas. While, if one imagines the very huge 
plasma gas corotating in the exterior and takes the limit for- 
mally L — > 00, then the wave number of region (3) should be 
equal to that of region (2), fc' 3 ' — ► . We thus have that 
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x+ — > 2k± , x_ — ► and 7 — > 1, which yield the boundary 
condition d/dz[\n^^] = ik± . 

One can draw some important facts from the boundary 
condition 13511 . Especially when the wave number satisfies 
the periodic condition x±L — nn (n — 1, 2, • • • ), this yields 
7 = k^ /k± and thereby the boundary condition for pure 

vacuum is recovered d/dz[ln£±^] — ik^ . The background 
rotation has no influence only on the Alfven waves satisfied 
with this condition. Whereas, when the wave number satis- 
fies the another periodic condition x±L = (2n— l)i"/2 (n = 
1, 2, ■ ■ • ), one obtains 7 = k± /k^ and thus the boundary 
condition works out to be d/dz[ln^ } ] = ik (3) /fc (3) ) 2 > 



ik^ 3 \ The RHS of this expression turns out that the waves 
satisfied with this condition are highly transmitted, as far 
as fc±' does not change much over one wavelength, say the 

validity of WKB approximation fl/&£^l jd/c^'/dzj <C fc£^ 
holds. This peculiar result means that the internal shear 
modes can be coupled onto the high frequency Alfven modes 
at the stellar surface, which drastically enhances their trans- 
mission rate at certain wave numbers or frequencies. 
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